Photon number discrimination without a photon counter and its apphcation to 

reconstructing non-Gaussian states 
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The non-hnearity of a conditional photon-counting measurement can be used to 'de-Gaussify' a 
Gaussian state of light. Here we present and experimentally demonstrate a technique for photon 
number resolution using only homodyne detection. We then apply this technique to inform a 
conditional measurement; unambiguously reconstructing the statistics of the non-Gaussian one and 
two photon subtracted squeezed vacuum states. Although our photon number measurement relies 
on ensemble averages and cannot be used to prepare non-Gaussian states of light, its high efficiency, 
photon number resolving capabilities, and compatibility with the telecommunications band make it 
suitable for quantum information tasks relying on the outcomes of mean values. 



Historically, quantum optics has been divided into two 
complementary camps, each exploiting only one aspect 
of the wave-particle duality of light. The "Discrete Vari- 
able" (DV) approach relics on quantized measurements 
of optical systems with low photon numbers, where the 
qubit space is often spanned by two orthogonal polar- 
izations. In contrast, the "Continuous Variable" (CV) 
approach focuses on field measurements of compara- 
tively bright beams, with states defined in an infinite- 
dimensional Hilbert space. With quantum optics provid- 
ing an ideal field for the first demonstrations of quantum 
information protocols, non-Gaussian states and transfor- 
mations are a necessity, forming a crucial resource for 
quantum communications [l|, 0], metrology 01 and 
quantum computing 0-0] ■ techniques are popular 
candidate quantum information protocols, owing to their 
high detection efficiency and compatibility with existing 
telecommunication infrastructures. However, the usual 
continuous variable toolbox, comprising Gaussian states, 
linear optics and homodyne detection, is insufficient to 
break out of the Gaussian regime. 

In the absence of extreme nonlinearities, integrat- 
ing DV photon-counting into continuous variable se- 
tups has been shown to provide probabilistic 'de- 
Gaussification' [8l-[l5|. This approach circumvents the 
'Gaussian' limitations of CV whilst enjoying the degrees 
of freedom both techniques provide. However, these 'hy- 
brid' experiments face challenges arising from simultane- 
ously exploiting both the wave and particle properties of 
light. In this paper we demonstrate a continuous vari- 
able analog to the 'photon counter' that enables discrim- 
ination of quanta of light. We show the versatility and 
efficacy of this technique by experimentally reconstruct- 
ing the non-Gaussian one and two photon subtracted 
squeezed vacuum (1 and 2-PSSV) states. 

CV techniques combined with linear optics are known 
to be insufficient to prepare non-Gaussian states with 
negativity in their Wigner functions [l^. Nevertheless, 



the idea of measuring the corpuscular nature of light with 
only CV techriiques has been theoretically [l9l - l2]| and ex- 
perimentally j23 - l24| investigated. Here we extend these 
ideas and show how DV heralding can be replaced by pure 
CV conditioning for the reconstruction of non-Gaussian 
states. This protocol avoids experimental issues arising 
from 'hybridising' a setup [l^ whilst harnessing the ex- 
isting mode selectivity, high quantum efficiency and low 
dark noise of homodyne detection. Although the fc-PSSV 
states are not heralded, remarkably we can still extract 
their quantum statistics. Using this method, we have 
successfully reconstructed the 1 and 2 PSSV states. 

A beam splitter is used to divert a small portion of an 
input squeezed vacuum (mode a) for conditioning whilst 
the remainder (mode b) is sent to a tomographic homo- 
dyne detection that measures = e^*^at + e*^aj, where 
cib and al are the annihilation and creation operators in 
mode b and 9 is the quadrature angle. The probability 
distributions of the fc-PSSV are given by 



4{x) = 7r{Xli=x\na = k). 



(1) 



In hybrid experiments (see FigUJi)) these are simply es- 
timated by reconstructing histograms of only when k 
photons are detected in mode a. The CV-only approach 
presented here (see FigHJii)) replaces the DV condition- 
ing with a dual-homodync measurement and exploits the 
relationship between photon number and field operators: 



a^aa = 



{xtf + {x-J^ 



(2) 



where X+ = A^f ~° and X^ = xt '^^^ , or any other pair of 

orthogonal quadratures A^ and xt^'"^'^ ■ Perfect simul- 
taneous measurement of conjugate observables A+ and 
A^ being forbidden, the two observables A^^ and A^ 
sampled by the dual homodyne detection incur a noise 
penalty in the form of vacuum fiuctuations entering the 
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FIG. 1: (Color online) Schematic of fc-photon sub- 
tracted squeezed vacuum states reconstruction. Both 
schemes require squeezed Gaussian input states. A beam- 
splitter is used to tap-off a small portion of the input beam, a, 
for conditioning measurement. The remaining light, b, is anal- 
ysed with tomographic measurements, (i) Hybrid Setup: 
outcomes from photon number resolving detectors are used to 
gate the results of the tomographic detection. Keeping only 
the conditionally heralded statistics, non-Gaussian quantum 
states are reconstructed; (ii) Pure Continuous Variable 
Setup: simultaneous orthogonal quadrature measurements 
are performed to weight the statistics of the tomographic mea- 
surements, producing non-Gaussian statistics. 



empty input port v of the beam-splitter, (see FiglTJii)): 



Xti = {X++X+)/V2 
X-^ = {X--X-)/V2. 



(3) 



As a consequence, attempts to measure Ua by using 
dual homodyne detection produce continuous real val- 
ues rather than integer results. This prevents any direct 
heralding of the fc-PSSV with Gaussian measurements 
only, in accordance with [Tsj , as exploiting the measure- 
ment outcome for single shot heralding would lead to 
meaningless results. However, since a tomographic recon- 
struction only deals with ensemble measurements, and 
because the statistical properties of the vacuum fluctu- 
ations are perfectly known, this additional noise can be 
effectively cancelled. We first rewrite 7r^(a;) as 



TT,, X 



7r(X« = a-, = k) (5(X« - x)5i,^k) 



n{na = k) 



7r{na = k) 



(4) 



We then note that Jft^j, can be approximated by a poly- 
nomial Pkifia) equal to zero at any significant photon 
number other than fc, for example Pk{na) — f^ai^a — 
l)...{fia — k + 1) in the limit of low squeezing. Dropping 
the normalisation factor, we find 



T,i{x),^{5{X',~x)Pk{ha)). 



(5) 



This means 7r^(a;) can be experimentally obtained by con- 
structing histograms of X^ where all values are taken 




FIG. 2: (Color online) Experimental set-up for the pure 

CV protocol. A squeezed vacuum state (Vs = —3.78 dB of 
squeezing for Va ~ -1-4.33 dB of antisqueezing) at the side- 
band frequencies between 3.6-4.2 MHz of a bright carrier is 
produced by a degenerate optical parametric amplifier (OPA). 
A variable beam splitter, a A/2 wave plate combined with a 
polarising beam splitter (PBS), reflects 20% of the incom- 
ing beam for conditioning and transmits the remaining for 
tomographic reconstruction. The conditioning beam is split 
between two homodyne detections (channels ai and 02) mea- 
suring two arbitrary orthogonal quadratures X^j and X^2 ' ■ 
The orthogonality between the homodyning angles for modes 
al and a2 is set using a combination of polarisation optics, 
whilst the global phase (p does not require active control (see 
eqn[2]). The tomographic reconstruction is preformed by sam- 
pling at 12 fixed angles 6 between and 180° in intervals 
of 15° . The experiment is controlled and automated using an 
FPGA based system 



but weighted by the corresponding value of P{na), ef- 
fectively discarding results corresponding to unwanted 
Ua ^ k. In view to reconstructing the statistics of 
the 1-PSSV state, the simplest applied weighting, where 
P[na) = iT-a can be understood as eliminating contribu- 
tions of X^ corresponding to = {i.e. no photon sub- 
traction), and retaining the desired contributions where 
Ua = 1 {i.e. successful photon subtraction). This triv- 
ial remark becomes powerful when injecting ([2) in ([S]), 
thus defining P{X'^ ,X~): wc can now use the real val- 
ued outcomes of our dual homodyne detection to weigh 
the results of our tomographic measurements. To do so 
we need to transform P{X'^ ,X^) into another polyno- 
mial Q{X^T^ X~2) of the measured observables that per- 
forms, on average, the same weighting. Q(X^j^, X~2), 
in general, is not trivial to obtain due to non commut- 
ing algebra. The demonstration of the existence of Q 
for any P{iia), and algorithms to compute it, will be 
published elsewhere. Here, as an illustration of the con- 
cepts used in the general case, we only prove the result 
Q{Xt,,X-2) = [(X+)2 + (l-)2 - 2]/2 for the simpler 
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(b) ha{ha - 2){ha - 3) 




(d) haiha - ^){n^ - 3) 




FIG. 3: (Color online) Wigner Functions reconstructed via 
the inverse Radon transform. The specific weighting polyno- 
mial is given at the top of each diagram. 1-PSSV: (a) with- 
out and (b) with higher order correction respectively, showing 
a improvement in the negativity from —0.018 to —0.055 com- 
pared to —0.16 in the ideal case. 2-PSSV: (c) without and 
(d) with higher order correction respectively. It shows an im- 
provement in the size of the central positive fringes (and first 
negative fringe) from -1-0.077 (-0.011) to -1-0.122 (-0.023) 
compared to -1-0.16 (—0.03) in the ideal case. The insets dis- 
play the corresponding calculated 1 and 2-PSSV states assum- 
ine nure initial soue.ezed states and ideal nhoton subtraction. 



case Pifia) = fia- Using ([2]) and ([3]) wc find 

(6{Xl-x)Q{X+„X-2)) = 

(5{Xt-x)n}j + {5{Xl-x)K) 

where 

A = (A+ + A-)/2 (7) 
A± = {X^f-l±Xtxt±XtX^. (8) 

We now note that the vacuum fluctuations X^ are un- 
correlated to those of X'^ and X^ and lience can be 
averaged out separately in ([6]), for example {6{Xjj — 
x)X+X+) = {S{X^ - x)X+){X+). Finally, injecting the 
values (Xjf) = and ((X^)^) = 1 given by quantum 
theory, we easily find that {S{X^ — 2;)A) = 0. Hence we 
can experimentally obtain tt^{x) through weighted his- 
tograms: 

= {S{X', - x)Q{X+,X-,)) . (9) 

As hinted at earlier, this idea can be extended to more 
sophisticated polynomials P{na) that approach photon 
number selectivity, permitting us to reconstruct purer 
and/or larger fc-PSSV states. Larger fc-PSSV states 
can be reconstructed using polynomials that remove the 
contributions associated with subtracted photon num- 
bers smaller than k, and keep outcomes corresponding 
io Ua = k as the predominant contribution. For ex- 
ample, the 2-PSSV state can be reconstructed using 
Pifi-a) = fi-aifi-a ^ 1), removing contributions resulting 
from having or 1 photon in mode a. 

Wc can also apply the same technique to remove contri- 
butions corresponding to unwanted higher order photon 
number subtractions contaminating the weighted ensem- 
ble. If we consider the simplest P{na) = na condition- 
ing polynomial to reconstruct the 1-PSSV, whilst con- 
tributions corresponding to ria = are cancelled, mea- 
surements associated with two photon subtraction are 
kept, and their statistical contribution is weighted at 2. 
As a result the reconstructed state is a statistical mix- 
ture of the 1-PSSV with some contribution from the 2- 
PSSV, producing a partial wash-out of the negativity 
of the reconstructed Wigner function. However, using 
P{na) = naifia — 2) for conditioning allows us to re- 
move contributions resulting from having or 2 photons 
in mode a. These two conditioning techniques can, in 
theory, be extended to an arbitrary order allowing a CV 
analog to the DV photon number resolving detector. 

Our experimental setup is detailed in Figl2] All of 
the presented Wigner functions (FigI3|) are reconstructed 
directly from the probability distributions obtained by 
applying various conditioning polynomials to a unique 
dataset. This dataset is composed of approximately 1.2 
billion samples for each detection mode and for each of 
the 12 tomographic angles. Whilst some reconstruction 
methods require assumptions on the nature of the state. 
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the inverse Radon transform [17[ used here is direct and 
assumption free. Moreover we do not correct for any 
experimental inefficiencies. 

FiglSUa) shows the Wigner function obtained using the 
conditioning polynomial P{na) = "-o, and displays a clear 
negativity. Our protocol, however, relies on correlations 
between modes al, a2 and 6, which are affected by any 
process that adds extra uncorrclated classical or quantum 
fluctuations. The primary limitation is the finite purity 
of our squeezed vacuum ressource, that can be modelled 
as an effective loss of ~12% applied on a pure squeezed 
state. We note that this non-ideal resource is also a lim- 
iting factor in hybrid experiments [l^ . Another possible 
limitation comes from the finite homodyne efficiency and 
dark noise. The dark noise, analogous to the dark counts 
of DV detectors, is at least 22 dB smaller than the vac- 
uum fluctuations. Our homodyne efficiency is limited 
by the quantum efficiency of the detectors, estimated at 
between 93-96%, whilst mode matching efficiency is typ- 
ically greater than 99%. In contrast, hybrid experiments 
suffer a substantial loss contribution from the difficulty 
in isolating the correct spatio-temporal mode in the con- 
ditioning and characterisation stages. Finally, the need 
for a finite conditioning tap-off and squeezing inevitably 
introduces spurious higher order photon subtraction con- 
tributions in the reconstructed state. By using the con- 
ditioning polynomial P{na) = fia{na — 2)(na — 3), the 
contributions from 2 and 3 photon subtractions are re- 
moved; the further higher orders being negligible. As 
expected, Figl3l[b) shows a considerable improvement in 
the negativity of the reconstructed 1-PSSV state. 

Figl3i;c) shows the 2-PSSV obtained using P{na) = 
"■a('«a — !)■ The reconstructed state exhibits the expected 
central positive fringe and two negative side fringes, as 



well as a bigger separation of the two coherent compo- 
nents [l6| . As shown in Fig|3|d), correcting for higher or- 
der contaminations by using P{na) — fiaifia ~ ^){na — 3) 
enhances the purity of the 2-PSSV state, evidenced by 
an improvement in the size of the fringes. 

In this paper we have demonstrated a photon dis- 
crimination technique based on a dual homodyne detec- 
tion, and used it to reconstruct conditional non-Gaussian 
states using only Gaussian resources and measurements 
and linear optics. Whilst this technique does not allow 
us to directly prepare non-Gaussian states, it enables 
their characterisation, and can be extended to charac- 
terise other quantum information protocols relying on 
hybrid, or 'de-Gaussification', techniques. By avoiding 
direct photon-counting, we circumvented the difficulties 
arising from simultaneously exploiting the wave and par- 
ticle nature of light, permitting us to unambiguously re- 
construct the I and 2-PSSV states. This idea of probing 
the 'quantised' nature of the quantum system via prob- 
ing of its continuous variables could also prove interesting 
for fields such as opto-mechanics, where direct measure- 
ments of the quantisation are unavailable or technically 
difficult. 
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